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A five- dimensional lattice space is decomposed into a number of four-dimensional lattices called 
'layers'. We can interpret the five-dimensional gauge theory as four-dimensional gauge theories on the 
multi-layer with interactions between neighboring layers. In the theory, there exist two independent 
, coupling constants: fli controls the dynamics inside a layer, and /3s does the strength of the inter- 

layer interaction. We propose the new possibility to realize the continuum limit of a five-dimensional 
theory using four-dimensional dynamics with large j3± and small /3s. Our result is also related to the 
03 higher dimensional theory by deconstruction approach. 
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1 Introduction 



Kaluza and Klein (KK)'s idea was a nice unification scenario that a four-dimensional gravity and a 



Maxwell theory are derived from a five-dimensional gravity [1]. But their approach was done within 
a classical framework and it is not so easy to realize its quantum theory. The first problem is that 
five-dimensional theories are generally unrenormalizable even if they are renormalizable theories in four 
dimensions. The another point is that the property of phase transitions in five-dimensional lattice theories 
is exactly same as that of mean-field theories. The latter implies that it is generally difficult to take the 
continuum limit in the five-dimensional theory using the second order phase transition and the associated 
fixed point (FP). Actually, some searches for both properties in a five-dimensional gauge theory have been 
tried [2, 3, 4, 5] and it is still far that one realizes a continuum limit from five-dimensional statistical 
mechanics. 

Although most of possibilities in the construction of five-dimensional theories are shut, there is a 
path to the construction, that is to use well-known properties of a four-dimensional statistical mechanics. 
Since we know the continuum limit of a four-dimensional gauge theory, the limiting operation keeping 
fixed physical quantities can be considered in the meaning of a four-dimensional theory. Recently, a new 
approach for the KK's scenario called deconstruction has appeared using four-dimensional asymptotically 
free gauge theories in the continuum [6, 7]. The approach realizes KK-like modes although they are 
completely four-dimensional theories. It is an interesting idea but can not give any answer for difficulties 
of higher dimensional theories because the idea is based on the perturbative picture of four-dimensional 
theories and has nothing to any clue about properties of the statistical mechanics. 
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On the other hand, Fu and Nielsen suggested to reduce the extra-dimensional space by the dynamics 
of a higher-dimensional lattice gauge theory[8]. In the idea, they considered a decomposition of a five- 
dimensional lattice space into an ensemble of four-dimensional lattice spaces (multi-layer) . Unfortunately, 
they studied about compact U(l) gauge theories which have no continuum limit in four dimensions[9, 
10, 11, 12]. In this paper, we formulate a five-dimensional lattice SU{2) gauge theory as a block of four- 
dimcnsional theories with an extra interaction like as Fu and Nielsen. The phase structure of the theory 
is clarified with Monte Carlo simulation and we discuss about the continuum limit as a four-dimensional 
multi-layer system (multi-layer world). Finally, we suggest how to construct a five-dimensional theory by 
explicit limiting operations. A limiting curve corresponding to five-dimensional scale shall becomes to an 
envelope. 

The outline of our paper is as follows: our model setting and the notation are explained in the next 
section. In sec. 3, our model is investigated by a numerical method. In sec. 4, we discuss physical properties 
of our model such as gauge squared mass matrix, the four-dimensional continuum limit and suggestion 
of a new idea for defining a five-dimensional theory. Sec. 5 is devoted to summary and discussions. 



2 Our Model 



In this section, we shall set our model in addition to explanation for a layer structure in the five- 
dimensional lattice. Link variables along a fifth direction are interpreted as scalar fields in the model. 
Naively, this has a picture of a four-dimensional theory in the continuum limit. Afterwards, measured 
order parameters in our numerical simulation are listed and an expected phase diagram is sketched. 

1. Layer structure and Wilson action Wc consider the five-dimensional pure 5/7(2) Wilson action 
on a Nf x N$ lattice with a periodic boundary condition: 



Slat = 04 ^ 



1-2-TrL^ 



(1) 



Trt/f 



= TV /7™ 5 77™ 5 rrt«s 7rt™5 
— 11 ^n,n u n+fi,v u n+i>,fj, u n,v 



= TrTT n ~> TT n ~> 7vt™5 + l 
— u n,n u n+[t,5 u n,n 



(2) 
(3) 



where the sum is defined as Ep 4 = E„,„ B £i<,.<„<4 » Ep 5 = E„,„ 5 £i<„<4 ■ According to Fu 
and Nielsen [8], we call the four-dimensional subspace a 'layer', because the individual four-dimensional 
subspacc must be independent of each other if (3§ is equal to zero. Every layer was aligned along the fifth 
direction — see Fig.l. An 715-th layer is noted £„ 5 . ?7™ 5 is a link variable inside £„ B and t/™ 5 5 connects 
between £„ 5 and C ns +i- These link variables behave as 2 of SU(2). Then the first term of the action (1) 
represents interactions inside layers controlled by fii, and the second term represents interactions between 
different layers (inter-layer) controlled by (3 5 * . Note that for convenience two coupling constants, 



/3= yfWs, J=\lf 

P4 



(4) 



are also used. 

2. Measuring quantities To understand the dynamics of (1), the following quantities are measured. 

• Expectation values of the plaquettc Up 4 inside a layer and U p 5 between neighboring layers. 

• Creutz ratio of loops inside a layer, 



1,Z„-1) ] 



W(l^l v -l)W(l 

where W(l^, l v ) implies an expectation value of Wilson loop with Z M x l v size on the same layer. 



(5) 



*A similar model with the adjoint Higgs field is invcstigatcd[13], where the field plays an essential role in the existence 
of a 'layered' phase. Although we does not assume the phase in this work, it seems an interesting subject. 
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• Polyakov loop along fifth direction, 



L 5 ,„^i T rn^; 5 5) . (6) 



2 

"5 

• Inverse of correlation length for Polyakov loops, 

(L 5 ,nL5, m ) - (L 5 , n ) 2 ~ cosh(-i? \n - m|) . (7) 

Up 4 and Up s correspond to parts of the internal energy. 04 (7 M , l v ) is an order parameter of four-dimensional 
confinement inside a layer and (£5,™) is that of Z 2 center symmetry t. 

3. Intuitive understanding of dynamics Depending on the magnitude of (34,05, the dynamics (1) 
has various feature. In the case that 7 is larger than unity, inter-layer interactions are strong compared 
to those inside layers. In this region the influence among the layer can not be ignored because a layer is 
strongly bound to its neighboring layers. It seems difficult to control the system in taking the continuum 
limit because there remain the strong extra-interactions caused by the fifth direction in the coupling 
region. As 7 is smaller than unity, interactions of inter-layer are weaker than those of inside layers. Then 
layers are loosely bound to each other, and in this region the dynamics is similar to a four-dimensional 
gauge theory with extra fields. We call this picture a 'multi-layer world'. 

From two order parameters 04 and (L 5 ) we can imagine four phases. Since our model is reduced to a 
four-dimensional gauge theory at /3s = 0, we know that the phase of cr 4 = and (L5) = is not realized 
— see Fig. 2. 

4. Naive continuum limit as four-dimensional gauge theory We consider the naive continuum 
limit of our model for a particular case which ?7™ 5 5 is interpreted as a fundamental 'scalar' field V™ 5 on 
the Cn 5 ■ In the interpretation, the fifth direction is simply an internal degree of freedom. We should 
consider limit CI4 — > as the naive continuum limit because there is only lattice spacing 04 inside layers. 
In the first term of the action (1) S„ 5 corresponds to a trace operation in the internal space and S ra dose 
to the summation of four-dimensional space. Then the relation between the inside layer coupling (3^ and 
a four-dimensional gauge coupling 34 is determined as 

& = 4 ■ (8) 

On the other hand, the second term of the action (1) means interactions between gauge field and 'scalar' 
field written as 

Up 5 = U^V^.u^+iy^. (9) 

The /?5 is the inter-layer coupling which controls dynamics of 'scalar' fields, and it is determined from 
the renormarization factor in the continuum limit. 



3 Numerical Analysis — Phase Structures on Various N 5 — 

In this section, we investigate the phase structure of our model by using Monte Carlo simulation. We 
simulate for Wilson action (1) using Nf x N 5 lattice (iV 4 = 8,10,12, N 5 = 1,2,3,5,6,8) by heat-bath 
algorithm. We measure expectation values of plaquettes (Up 4 ), (£7p 5 ), specific heats C4 = (f/p ) — (Up 4 ) 2 , 

C5 = (Up ) — (f/p 5 ) 2 , Creutz ratio and Polyakov loop with its fluctuation X5 = (^5) ~ (^s) 2 - The 
statistical errors of these measurements are estimated by the jackknife method. In the following figures, 
open and filled symbols are noted as results of ordered and disordered start respectively. We describe 
N 5 — 1 and N^> \ cases separately. 

tFor the JV5 = 1 case, a singlet part of U n ,5 is important for the phase structure. For the realization of a non-trivial 
transformation for the singlet part only, we add a special gauge transformation to the center symmetry. 



3 



3.1 Case of iV 5 = 1 



Fig. 3 represents the phase diagram at N$ = 1 and plots measured lines where (i) 04 = 0, (ii) 7 = l/\/l2 <~ 
0.29, (iii) fa = 10.0 and (iv) 5 = 10.0. Our calculations support the expected phase diagram (Fig. 2), 
and their detailed discussions are summarized as follows: 

(i) 04 = 0: Fig. 4 shows results of (L 5 ) and X5- Those quantities imply that the confincment- 
deconfinement phase transition occurs at 0^ c <~ 0.6. Especially in x 5 , we can see the volume effect 
which is a characteristic signal of the phase transition. 

(ii) 0^ = 10.0: From (L 5 ), the similar signal of the transition can be seen at 0$ ~ 0.6. In addition, 
(T4 also bends sharply at 0$ <~ 0.6 — see Fig. 5. 

(iii) 05 = 0.3: In Fig. 6, we can confirm that signal of crossover about 04 is exactly same as that of 
an ordinary four-dimensional lattice gauge theory (/3 5 = 0). 

(iv) 05 — 10.0: We can see that 04 is decreasing as 04 is growing in Fig. 7, although it is difficult to 
recognize a rapid change of 174 at 04 ~ 1.2. 

The naive continuum limit of an N 5 = 1 system corresponds to the four-dimensional gauged non-linear 
a model: 

Slat- J d i x^Fl v+ l -fl(D,V{x)) { (D,V{x))^ (10) 

where V(x) is a scalar field with a constraint condition 

V(x) = 6 {x) + ia ■ 6{x) , 

el + e 2 = \. (11) 

This condition reads to (<9 M (9o) 2 <C (d fl 9 a ) 2 at 9 a <C 1. Then the system corresponds to a Georgi-Glashow 
model with a Higgs triplet [14, 15, 16] which is equivalent to a four-dimensional 0(3) Heisenberg model. 

3.2 Case of iV 5 > 1 

The phase diagram for N 5 > 1 and measuring directions are shown in Fig. 8. Dashed lines are (i) 7 = \/2, 
(ii) 7 = 0.5, 1/V2 (~ 0.71), 1/1.2 (~ 0.83), (iii) 7 = 0.29 and (iv) 4 = 8.0, 10.0 and 20.0 respectively. 

(i) 7 = \f~2~: In Fig. 9, as becomes to large, the new phase appears at <~ 1.2 where the phase varies 
from (74 7^ and (L5) = to (74 7^ and (L 5 ) 7^ 0. Afterward the new phase transition occurs at <~ 2.0 
where the phase varies from (74 7^ and (L 5 ) 7^ to (74 = and (L 5 ) 7^ 0. 

(ii-a) 0.5 < 7 < 0.83: In the middle range of 7, we found a first order phase transition for both order 
parameters Up 4 and Up 5 on a 8 4 x 2 lattice in Fig. 10. There the hysteresis loop is clearly seen at 7 = 0.71. 
As N 5 approaches to N4, the range is wider and includes unity where the theory is recovered as usual 
five-dimensional one. Then this signal of the first order transition corresponds to that of the isotropic 
five-dimensional SU(2) gauge theory [3]. 

(ii-b) 0.0 < 05 < 1.2: Fig. 11 shows a change of (74(2, 2) for 04. To compare with it, the result of an 
ordinary four-dimensional lattice is also shown at the circle symbol. As 05 is larger, a slope of (74(2, 2) is 
deeper. The signal of crossover is changed to the second order phase transition at 05 ~ 1.0 and further 
the first order phase transition at 5 ~ 1.2. The change of (74(2, 2) for 5 at 04 = 2.3, is shown in Fig. 12. 
We can see that (74(2, 2) or the tension is slightly varied in small 5 , although it becomes clear in 5 > 0.8 
near the phase transition. 

(iii) 7 = 0.29: We can find the phase transition at ~ 2.3 (04 ~ 0.79, 0$ ~ 0.67). It is a surprising 
fact in a multi-layer world that the transition point is independent of N 5 > 1. Fig. 13 shows results of Up 5 
and C5 on a 8 4 x 7V 5 lattice, where N 5 = 1, 2, 3, 5, 6 and 8. Both quantities for N 5 > 1 are unchanged. 
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(iv) = 8.0, 10.0 and 20.0: At /3 4 = 10.0 we can see that the phase transition occurs at /3 5 <~ 0.63 
from (74 7^ and (L 5 ) = to 04 = and (L 5 ) ^ — see Fig. 14. The hysteresis loop can not be seen 
in their signals. Therefore we expect the phase transition is second order. The inverse correlation length 
for Polyakov loops shows at /3 4 = 10.0 on a 12 4 x 3 lattice in Fig. 15 and also supports the property of the 
phase transition. In addition, we found that the transition point for N 5 > 1 is closer to that for N 5 = 1 
as P4 is larger — see Fig. 16. 

Those interesting phenomenon in the multi-layer world shall support our discussions in the next 
section. 



4 Multi-Layer World 



4.1 Multi-layer world and gauge fields 

In this section, we investigate a multi-layer world (7 is sufficiently small to unity and iV 5 > 1). This 
world is composed of four-dimensional gauge theories bound by 'scalar' fields: = 5 . The fields are 
transformed as bi- fundamental representations of and A+i- The gauge transformations are defined as 



Y l (n) Y % (n + (1) for a gauge field inside a £j, 
V l n -» Y l {n) V l n Y l+1 \n) for a 'scalar' field between d and £ i+1 , 

where Y l {n) is an arbitrary element of SU(2). From (13), V£ can be fixed to unity except for V„, 



(12) 
(13) 



V 1 ' 

v 2 ' 



V K 



= Y^n) V* Y 2 \n) 
= Y 2 {n) V 2 Y 3 \n) 



= 1 
= 1 



(14) 



= Y N *(n) V™* Y l \n) = Y l (n) \\[^ F lf (n) = Y\n)X n Y 1 \n) . 



Note that TrV^ 5 is exactly a Polyakov loop L§ n . Since a non-trivial can be moved to every inter- 
layer, this implies that the system has a translational invariance to the fifth direction. In this point, the 
multi-layer world is significantly different from a just four-dimensional theory. 

Based on this setting of the gauge fixing, we consider mass spectrum of vector fields, A® % . The masses 
include the effect of fluctuations by layers. We get the following mass term, 



N 5 



■ • 1 2 ' 

n,/i i,j=l 



(15) 



where the mass squared matrix is defined as 

1 -1 



M ■ 



4 a 2 



2 

-1 



L 



5,n 



-1 



4a 2 , lJ 



(16) 



using Un p = exp {i a 4 34 ^- A^ ^} . For the purpose that we go ahead to non-perturbative analysis, the 
mass may be also non-perturbatively defined by the correlation length between different layers. But the 
analysis is necessary for a large N 5 lattice and it is beyond our present computer capacity. Therefore, 
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we take the semi-classical approach such that the mass formula is derived by expanding a field A n ^ and 
the mass squared matrix is calculated by Monte Carlo simulation. Under this approximation, L 5 , n can 
be replaced with (L5) = 1/(214) ■ XL Tr A„, an d V4 is the four-dimensional volume of a layer. In the 
case of (£ 5 ) 2 = 1, we can easily analyze this matrix Mij owing to the periodic boundary condition for 
the fifth direction [6, 7]. Eigenvalues of this matrix are 



2 



A fe = 4sin^— , fc = 0,l,---,JV 5 -l , (17) 
^ * 5 



and at k -C N5, they are simply written as 



2irk 



2 



A fe =(^H • (18) 



It is noted that the values for even k are independent of (L5). Since it is not so easy to get the 
values for odd k, we numerically calculate them with arbitrary (L5) which is determined by Monte Carlo 
simulation. Fig. 17 shows plots of eigenvalues A^ versus (L 5 ) 2 at N 5 — 3, 6 and 12. 

We see that eigenvalues are doubly degenerate at (L5) 2 = 1. In addition, a noteworthy thing is that 
the eigenvalues at (L 5 ) = 1 and N 5 = n are equal to the eigenvalues at (L 5 ) = and N 5 = n/2 except for 
its degeneracy. Our system at (L 5 ) 2 = 1 and N 5 — n has really the corresponding relation to a circular 
model at N5 = n, and a system with {L5) = and N5 = n/2 has a relation to a model with Z2 orbifolding 
[6, 7]. This circumstance means that S 1 and S 1 /Z2 with respectively different topologies are connected 
by one parameter (i 5 ) which corresponds to the energy releasing quarks from a layer. The reason is the 
existence of 'singular' gauge transformations such as (14) on lattice theories in contrast to continuum 
theories. The eigenvalues for even k are topology-independent and mainly used in our following analysis 
because we do not have a clue to change to the topology although topology changing phenomena is an 
interesting topics. 



4.2 Continuum limit as a four-dimensional gauge theory 

In this subsection, we shall consider the continuum limit of four-dimensional gauge theories on layers. 
Since our fifth direction is identified with just a internal space, we take only the continuum limit 04 — > 
on layers. From /3 4 = 4/g|, we can express the mass squared matrix as 



M l = -^2^k({L 5 (P4,0 5 )),N 5 ) = ^A fe ((L 5 (/34,/3 5 )),^ 5 ) • (19) 
P4Q4 04 

At k <§C A^5, the simple mass formula 



is obtained. In order to take the continuum limit 04 — > with a fixed mass (19), a possible method exists; 

N 5 fixed, 7-^0 and — finite . 

04 

The approach is realized on /3 4 — > 00 with finite /?5. This theory is defined as four-dimensional asymptot- 
ically free one using 04 in the confinement phase like a la deconstruction. Actually, the scaling for 04 is 
non-trivial because we have to take care of bi-fundamental fields V// 5 . Our numerical result shows that 
the scaling is subtly different from an ordinary four-dimensional theory with (3^ — — see Fig.ll. 

Although wc can not use 04 in the case of the deconfinement phase, from (7) Higgs mass, 

H(0 4 ,p 5 ,N 5 ) 

m h = , (21) 

04 



G 



can be considered instead. This quantity is not yet seen as a good data in statistics in Fig. 15. But near 
f3 5c in both phases, it can be seen as the important scaling parameter 

H oc |/? 6 - /3 5c r (22) 

where v is a critical exponent for the correlation length and seems less than unity. It is consistent with 
second order phase transition at /3s = /?5 C ~ 0.6. The observation is supported by no hysteresis loop 
and the peak of the specific heat. In this standing point, I our -dimensional theories with inter-layer 
interactions are defined for every (a multi-layer world). 

The phase diagram for N5 > 1 in the multi-layer world is not changed by N$ — see Fig. 13. This fact 
suggests us an extension of a multi-layer theory. If we try to extend this internal space to an real extra 
dimension, N 5 must be took to infinity. When N 5 goes to infinity, vanishes from (20). To avoid it, 
it is necessary for a relation among 7,04 and 7V 5 . If we can introduce a new lattice spacing, 05 <~ 04/7, 
our squared mass (20) is written as 

This relation is exactly a KK mass. Originally, our fifth direction is just an internal space but (23) shows 
the possibility that the internal space is interpreted as fifth direction. Although the explicit definition 
for a 5 is discussed in the following subsection, we note that a relation 

implies just that anisotropy of gauge coupling constants corresponds to that of lattice spacings like as 
finite temperature theories. 



4.3 Suggestion for a new definition of five-dimensional gauge theory 

In the end of the previous subsection, we saw the possibility of an extra dimension. To construct a 
continuous extra dimension, we must consider the continuum limit (a 5 — > 0) combining with another 
limit 04 — > 0. This is a hard parameter-tuning but it becomes possible by the balance between extra- 
dimensional dynamics and four-dimensional one. In order to set the balance, we can define the quantity 
from (7), (19) and (22), 

ml (3 4 H*((3 4 ,(3 5 ,N 5 ) Nl(3 4 \f3 5 -f3 5c \ 2 » ■ [ °> 

If this dimensionless quantity is kept as a nonzero finite value in the limit, it is possible to make both a 
Higgs mass and KK one of the second excited mode finite. 

Let us consider when 1Z can be nonzero finite. From the a 4 — > and (25), we should investigate the 
region where 4 is large and H is small (/% ~ fee)- Our numerical study in section 3 suggests to us that 
the phase structure and the parameter /35 C unchanged under various N5 > 1. So, we can concentrate on 
large (3 4 and small H independent of N$. After of all, we can introduce a scale parameter 125, 

a5 55 v»7 4 - (26) 



From this parameter, we can define an extra-dimensional physical scale, 

R 5 = N 5 a 5 = -^a A . (27) 
VA27 

Can we realize this as a five-dimensional space? The answer shall be yes if the five-dimensional scale 

Rs = M^^ = ^~ (28) 
V P5 \/\imh vHrrih 
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is kept nonzero finite. From the relation (28), we find that R5 is larger as 04 is larger in the case of fixed 
05. The formed lines in Fig. 18 are expected to make an envelope in 05 ~ 5c , which the limiting line 
appears in the plot after 04 — > 00. The existence of the envelope shall be an important clue in the proof 
that a five-dimensional field theory is constructed with a finite extra-dimensional scale R5 and a finite 
Higgs mass. 

5 Summary and Discussions 

In this paper, we have investigated the phase diagram for a five-dimensional SU (2) gauge theory. From 
the view of a layer world a la Fu-Nielsen, this system is considered as an ensemble of four-dimensional 
gauge theories on layers with inter-layer interactions. In the region of small but nonzero 05, the phase 
structure including a second order phase transition is surprisingly unchanged with N 5 > 1. In the limiting 
operations, we consider a 4 — > firstly because we know the limit of the four-dimensional gauge theory 
well. After the construction of this multi-layer theory, we have sought the possibility that we can keep 
an extra-dimensional scale and a four-dimensional scale finite. 

From the universality argument, our theory on 04 — > 00 and 05 = /3 5c corresponds to a four- 
dimensional Ising spin model. The model has a second order phase transition with v > 0.5 [17, 18, 19, 20]. 
What happens for our model in the case of 05 C = ? In v > 0.5, we can obtain a finite extra-dimensional 
scale by the same procedure as above discussion. If v — 0.5, we need logarithmic corrections for non-trivial 
theories. 

There is another approach called deconstruction. It is the ensemble of asymptotic free theories and 
has KK-like modes in four dimensions. Based on the construction of a four-dimensional lattice gauge 
theory, we study the non-perturbative property such as the vacuum expectation value of a five-dimensional 
Polyakov loop. Both approaches seem complimentary to each other. 

As remained problems, after detailed numerical computations along to our scenario, it is natural to 
try (A) an extension to matter fields, (B) restoration of the rotational symmetry in the five dimensions. 
Although, in this paper, we have mainly investigated the small 05/04, we may try the connection to the 
large 05 104 to solve these problems. 
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ns-th layer 




Figure 1: Layer structure of five-dimensional lattice. 
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Figure 2: Intuitive understanding of phase diagram in (2.1). 
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Figure 3: The phase diagram for N$ = 1. Solid lines imply phase boundaries and a dash-dotted line 
is originated in the crossover of a four-dimensional gauge theory. Dashed lines correspond to measured 
directions. 
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Figure 4: (L 5 ) and \5 at At = 0. 
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Figure 5: (L 5 ) and a 4 (2,2) at (3 4 = 10.0. 
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Figure 6: 04 and (£5) at (3 5 = 0.3. 
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8*x1, p 5 =10.0 




Figure 7: <r 4 and (L 5 ) at /3 5 = 10.0. 
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Figure 8: The phase diagram for N 5 > 1. The meaning of lines is same as Fig. 3. 
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Figure 11: The circle symbol shows 0-4(2, 2) on an 8 4 lattice. The others are results on an 8 4 x 3 lattice. 
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Figure 12: er 4 (2,2) at f3 4 = 2.3. 
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Figure 14: (L 5 ) and er 4 at /3 4 = 10.0. 
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;ure 15: The lattice counterpart H for a Higgs mass at j3& — 10.0 on a 12 4 x 3 lattice. 
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Figure 16: U Ps at f3 4 = 8.0, 10.0 and 20.0. 



25 



N 5 =12 N 5 = 6 N 5 = 3 




Figure 17: Eigenvalues of a matrix Mij for N 5 = 12,6 and 3. 
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Figure 18: An expected envelope (fii = oo) formed by lines (3\, [3' 4 , (5'1, etc. 
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